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ABSTRACT: 

The effect of functionalization of carbon nanotubes on the thennal conductivity of 
nanocomposites has been studied using a multi-scale modeling approach. These results 
predict that grafting linear hydrocarbon chains to the surface of a single wall carbon 
nanotube with covalent chemical bonds should result in a significant increase in the 
thennal conductivity of these nanocomposites. This is due to the decrease in the 
interfacial thermal (Kapitza) resistance between the single wall carbon nanotube and the 
surrounding polymer matrix upon chemical functionalization. The nanocomposites 
studied here consist of single wall carbon nanotubes in a bulk poly(ethylene vinyl acetate) 
matrix. The nanotubes are functionalized by end-grafting linear hydrocarbon chains of 
varying length to the surface of the nanotube. The effect which this functionalization has 
on the interfacial thermal resistance is studied by molecular dynamics simulation. 
Interfacial thermal resistance values are calculated for a range of chemical grafting 
densities and with several chain lengths. These results are subsequently used in an 
analytical model to predict the resulting effect on the bulk thermal conductivity of the 
nanocomposite. 


1. Introduction 

Polymers have many applications and uses but are often limited by their inherent 
physical properties. Most polymers have very low thermal conductivity. If the thermal 
conductivity could be improved without sacrificing the other properties of the polymer, 
more applications would be available for these materials [1], A common commercial 
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polymer with many uses is poly(ethylene vinyl acetate) (EVA), which has a relatively 
low thermal conductivity (0.2 W/mK), typical for polymer materials. An approach of 
current interest to improve the thennal conductivity (TC) of these materials is the 
selective addition of nanostructured materials with high TC. 

Although many novel nanostructured materials are available for study, the 
material chosen for this study is a single-walled carbon nanotube (SWCNT) composite. 
Carbon nano tubes (CNT) are known to have extremely high TCs [2] and research in 
measurement and modeling of carbon nanotube composite materials is currently very 
active. Thermal conductivity has been measured for a number of CNT polymer 
composites. However, to date, only modest increases in total TC have been measured 
experimentally [3]. At the nano-scale, it is assumed that there are a variety of intrinsic 
material structures which may be influential in establishing the properties of these 
nanocomposite systems. Previous theoretical work [4-6] has suggested that the modesty 
of these TC increases could be attributed primarily to the presence of a large interfacial 
thennal (Kapitza) resistance between the CNT and the surrounding polymer matrix. It is 
assumed that this interfacial resistance can be directly influenced by the degree of 
functionalization at the interface of the CNT and the polymer matrix. Recent simulation 
work has indicated the influence of interfacial molecular structure on interfacial thermal 
resistance [7], 

The goal of the atomistic-based modeling presented herein is to address the issue 
of interfacial thermal resistance by exploring key structure-property relationships and 
establishing, through parametric studies, if increases in the TC of the composite can be 
achieved by proper functionalization of the CNT-polymer matrix interface. The 
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interfacial thermal resistance, R K , between the CNT and surrounding matrix is modeled 
with atomistic molecular dynamics (MD) simulations. Parameters associated with the 
functionalization of the CNT-EVA matrix interface are systematically varied. These 
parameters are the grafting density, cr, and length, n, of linear hydrocarbon chains 
covalently bonded to the CNT surface. The effects of defects in the CNT on the 
interfacial thermal resistance are also examined. Analytical modeling is also applied in 
order to consider the effect of various parameters on the TC of the composite. Key 
parameters to be varied in this analytical modeling include CNT length, CNT volume 
fraction, CNT thennal conductivity and interfacial thennal resistance. 

2, Method 

In the current work, computer simulations of SWCNTs in EVA have been 
performed. Atomistic structures with variations in CNT functionalization are first 
prepared using standard molecular dynamics (MD) techniques. A methodology [4-6] is 
applied to calculate the interfacial thermal resistance, Rk, in these structures. These Rk 
values are then used with analytical modeling in order to consider the effect of chemical 
functionalization of SWCNTs on the thermal properties of the resulting nanocomposite. 

2. 1 . Structure preparation 

In order to predict the Rk values between the SWCNT and the matrix, MD 
simulations were perfonned with short linear hydrocarbon chains (-C„H 2 „+i, n = 8, 12, 18, 
50) randomly end-grafted to the surface of a (6,6) SWCNT. The SWCNT was 
functionalized by this covalent end-grafting with a range of grafting densities (cr= 0-0. 1 1 
A' 2 ). Fig. 1 illustrates a typical SWCNT with end-grafted chains. The diameter of the 
SWCNT was approximately 1 nm. These functionalized SWCNTs were embedded in an 
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EVA matrix in the MD simulation. The amorphous EVA copolymer consisted of 15% 
vinyl acetate monomers by weight or 5.5% mole fraction. In these simulations, 70 chains 
of EVA copolymer containing 25 repeat units each were simulated as a polymer matrix 
surrounding a (6,6) SWCNT. The copolymer sequences for the EVA matrix were 
randomly constructed; however, this is relatively insignificant since each chain only 
contained 1-2 vinyl acetate monomers. The functionalized nanotubes were simulated in 
the bulk EVA matrix as a nanocomposite. Bulk density was effected through application 
of periodic boundary conditions 8 in a simulation box size of about 46 A x 46 A x 58 A. 
The SWCNT was extended through the z (axial) coordinate of length 58 A. In all the 
simulations, the PCFF force field 9 was used with the LAMMPS [10, 11] MD simulation 
software. The SWCNT atoms are not constrained in any way other than that necessitated 
by the periodic boundary conditions. 

A total of 17 structures were constructed and thermally equilibrated. These 
included one structure with no functionalization of the SWCNT and 13 structures with 
variable length, n, and grafting density, cr, of end-grafted chains. Three more structures 
were variations involving defects of the nanotube. The largest range of grafting densities 
was constructed with the 8 carbon chain (-C h H 2 «+i, n = 8) end-grafts, where the SWCNTs 
with the following four grafting densities were simulated in an EVA matrix: (cr= 0.0068, 
0.013, 0.054, 0.1 1 A' 2 ). The SWCNTs with 12 carbon chain end-grafts ( n = 12) were 
simulated at grafting densities: (cr= 0.0034, 0.0068, 0.013, 0.027, 0.054 A' 2 ). The 
SWCNTs with 18 carbon chain end-grafts (n =18) were simulated at grafting densities: 
(cr= 0.0068, 0.013 A' 2 ), and the SWCNTs with 50 carbon chain end-grafts (n = 50) were 
simulated at grafting densities: (cr= 0.0034, 0.0068 A' 2 ). 
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The bulk atomistic nanocomposite model is constructed by compression from low 
density. This procedure begins by placing the SWCNT with grafted chains and the 70 
polymer molecules in a box with periodic boundary conditions. This initial density is 
about 5% of the bulk density of 0.8 g/cm . To prevent the polymer chains from 
collapsing, the Van der Waals (VDW) potential was scaled by a factor of 0.001 and a 
short cutoff (3.0 A) was applied with Coulomb interactions turned off. To allow the 
polymer chains to diffuse through the periodic box and achieve random configurations, 
the MD simulations were run for 50 ps using a 1 fs timestep at 500 K, reducing the 
volume at a linear rate to a factor of about 10% of the final bulk density. A second MD 
run condensed the system to a density of about 50% of the final bulk density over 50 ps at 
500 K. The unsealed VDW potential is then applied with a cutoff of 9.0 A and Coulomb 
interactions are enacted with the Ewald summation technique. Energy minimization is 
performed on the system followed by constant pressure MD simulation. MD simulation 
is applied for 10 ps while increasing the pressure from 1 to 100 atms at 300 K. This is 
followed by a compression at 100 atms for 50 ps at 300 K. The pressure is reduced from 
100 to 1 atms over 10 ps. This is followed by a MD run at 1 atm for 100 ps and 300 K. 

2.2. Simulation of the heat transfer 

Following the last MD run of the atomistic structures at 300 K described in 
Section 2.1, a methodology [4-6] was applied to obtain the interfacial thennal resistance, 
Rk . The interfacial thennal resistance calculation involved the instantaneous heating of 
the SWCNT atoms to a temperature of 500 K. Following this instantaneous heating, the 
simulation was run at constant energy, as thennal energy was transfened from the 
SWCNT to the surrounding polymer matrix. The difference in temperature, AT, between 
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the SWCNT atoms and the other atoms was recorded as a function of MD timestep. Fig. 


2 shows a plot of the ln(A7) as a function of time since the application of the 
instantaneous heating. A least squares fit to this curve yields a slope which is the 
negative inverse of the characteristic decay time, z. The interfacial thennal resistance, Rk, 
[6] is calculated from Eq. 1, where c T /A T is the heat capacity per area of SWCNT surface 
(using a value of 5. 6x1 O' 4 J/nrK for cr/A r [4].) 


Rk = 


(c T / A t ) 


( 1 ) 


2.3. Analytical modeling of the thennal conductivity 

These interfacial resistance values from the heat transfer simulations are then used 
in an analytical fonnula [12] to calculate the TCs of nanocomposites. While molecular 
modeling yields Rk values, an analytical model is applied in order to calculate the thermal 
conductivity of the nanocomposite material. This model assumes a random orientation of 
the dispersed nanotubes within the polymer matrix and uses an effective medium 
approach. This approach yields a fonnula for calculating the TC of carbon nanotube 
composites [12], This is summarized by Eqs. 2-5. 
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Material properties K e , K c , and K m are the thermal conductivities of the nanocomposite, 
the CNT and the polymer matrix, respectively. In evaluating Eqs. 2-5, a value of 0.2 
W/mK is used for K m . Two values are used for K c : 6000 W/mK and 3000 W/mK [2, 5]. 
The diameter of the (6,6) SWCNT, d, is 8.136 A, while /is the volume fraction, L is the 
nanotube length and aK is the Kapitza radius. 

3. Results and Discussion 

This section will describe the results of the modeling and provide discussion on 
the observed behaviors. In Section 3.1 the prepared atomistic structures are examined. 

In Section 3.2 the results of the MD heat transfer are discussed. From these simulations, 
the interfacial thennal resistant values, Rk, for a range of functionalization conditions are 
calculated. In Section 3.3, these Rk values are then used with analytical modeling [12, 13] 
to predict the TC of the composites. By systematically varying the CNT length, volume 
fraction, thermal conductivity and the interfacial thermal resistance, key structure- 
property relationships were established and related directly to changes in the predicted 
TC of the composite. 

3.1. Radial density profiles of the prepared structures 

The radial density profiles of the nanocomposites are examined in order to 
determine if there are large differences in mass distribution, particularly near the CNT- 
matrix interface which may be significant. Fig. 3a shows the radial distribution functions 
of the mass density as a function of distance, r, about the center of the SWCNT. These 
profiles are shown for the MD simulations with n = 12 linear hydrocarbon end-grafted 
chains. The large peak at ~5 A is due to the CNT wall. At low to moderate grafting 
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densities, (cr= 0 - 0.027 A' 2 ), the density profiles show the familiar oscillating behavior 
associated with packing polymer chains at surfaces. At high grafting density (cr= 0.054 
A' 2 ), the profile has a qualitative change in shape, due to the exclusion of polymer matrix 
chains from the SWCNT surface. Fig. 3b shows the same profiles given only by the 
contribution from the SWCNT and the end-grafted chains. A qualitative change in the 
profile shape occurs at high grafting density (cr= 0.054 A' 2 ) due to the steric crowding of 
the grafted chains. 

3.2. MD heat transfer results 

The results obtained from the MD simulation described in Section 2.2 are 
presented here. The 17 separate structures representing the variation in grafting density 
and end-grafted chain length are each subjected to the methodology involving 
instantaneous heating of the SWCNT atoms. Decay plots similar to Fig. 2 are obtained 
for each structure and the interfacial resistance values, Rk, are calculated according to Eq. 
1. Fig. 4a shows a plot of Rk as a function of grafting density for the various end-grafted 
chain lengths (-C„H 2 „+i, n = 8, 12, 18, 50). Large decreases in R K (>80%) are observed 
with increasing grafting density for all chain lengths. A consistent correlation between 
chain length and decreasing R K is predicted for corresponding grafting densities. 

The grafting considered above assumes a pristine SWCNT with no defects and 
grafting taking place at sp3 carbons located adjacent to the same locations as the 
previously unfunctionalized sp2 carbon atoms, i.e. direct covalent bonding to a pristine 
nanotube. However, it is reasonably assumed that during synthesis bonding could occur 
at defect regions in the nanotube wall. To consider this type of structure, the bonding is 
altered in several cases to approximate the effect of defects adjacent to bonding sites. An 
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sp3 carbon is bonded to 4 atoms. In the case of the grafted chains, the sp3 carbon of the 
nanotube is bonded to the terminal carbon atom of the linear hydrocarbon and 3 other 
carbon atoms in the CNT. In order to approximate the effect of a defect, 2 of these sp3- 
sp2 carbon-carbon bonds are broken. This does not result in the removal of any atoms, 
although the bonding characteristics of atoms adjacent to the grafting point are 
substantially altered. This is illustrated in Fig. 4b, where a planar graphite section 
representing the CNT wall is shown. A hydrocarbon chain is bonded to the CNT wall. 

With these defects in place, the simulation is repeated for the n = 12 chains for 
several grafting densities. Results, presented in Fig. 4c, compare the interfacial thennal 
resistance of the pristine grafted nanotube with the resistance of the defect induced 
nanotube. Very little difference in interfacial thennal resistance is predicted between the 
two cases, despite the very different bonding characteristics in the vicinity of the grafted 
point. Specifically, in the defect induced nanotube, the polymer chain is essentially end- 
grafted to the edge of a hole in the nanotube. The results indicate that for the materials of 
this study, SWCNT defects will not have a major effect on the interfacial thennal 
resistance as influenced by the grafting of chains. 

The results of the heat transfer simulations are summarized in Table 1. The range 
of grafting densities, a, is indicated by the columns. The range of end-grafted chain 
lengths, n, is indicated by the rows. The interfacial thermal resistance values are placed 
in this anay in units of m“K/W (xlO ). The values in parentheses are for the simulations 
with the defects in the SWCNTs. Comparing grafting densities in terms of percentage of 
functionalized carbons, there are 576 carbon atoms in the SWCNT. By percentage of 
functionalized carbon atoms in the CNT, the grafting densities are: (0.86%, 1.74%, 
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3.47%, 6.9%, and 13.9%) corresponding to grafting densities (cr= 0.0034, 0.0068, 0.013, 
0.027,0.054, and 0.108 A' 2 .) 

3.3. Analytical modeling results 

Using the analytical modeling presented in Section 2.2 [12, 13], the thermal 
conductivity of the composite, K e (RiJ is plotted in Fig. 5a. Here, a value of 1x10 6 m (1 
micron) is used for the length, L and K c = 6000 W/mK. The volume fraction of 
nanotubes in the bulk material,/ ranges from 1-5% in these plots. This volume fraction 
is a parameter used in Eq. 2-5 and is not derived from any details of the MD simulations 
in Sections 3.1 and 3.2. The solid line is computed for a continuous range of values, 
while the symbols are marked for the specific Rk values obtained from the (n = 8) MD 
simulations. These results illustrate the strong influence that Rk has with this common 
nanotube length. Fig. 5b shows similar results plotted with a longer nanotube length ( L = 
10 microns). Fig. 5c shows results plotted with a much longer nanotube length (L = 100 
microns). Comparing Figs. 5a-c, it is clear that the Rk value is much more influential on 
TC at shorter CNT length. All three plots retain the same limiting values (Rk = 0), yet 
have strongly differing TC values at finite R K . Fig. 6 shows the influence of NT length 
more directly. The TC of the nanocomposite is plotted as a function of NT length for 
several values of Rk- The five finite values of Rk are taken from calculated values shown 
in Fig. 4a for n = 8, and the limiting value (Rk = 0) is included. The plots clearly indicate 
the diminishing effect of CNT length with decreasing Rk- 

Estimates for the thermal conductivity of CNTs vary considerably [2, 14-17]. The 
effect of functionalization may decrease the thermal conductivity [17]. We have 
arbitrarily used a fixed value of 6000 W/mK for the calculations presented in Figs. 5a-c 
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in order to focus on the effect of R K on the TC of the nanocomposite. Fig. 7 illustrates 
the influence of the thennal conductivity of the CNT ( K c ), by recalculating the values 
from Fig. 5c with a change of K c from 6000 W/mK to 3000 W/mK. This illustrates that a 
lower thermal conductivity of the nanotubes would still allow a significant improvement 
in the resulting TC of the nanocomposite. 

Fig. 8 shows the results from Fig. 5c as a function of grafting density, cr . This 
more clearly illustrates the dependency of the resulting TC of the nanocomposite on the 
grafting density of chains (n = 8) attached to the CNT. The dependence of the thermal 
conductivity on the chain length, n, is shown in Fig. 9. The grafting density for the data 
in this plot is cr = 0.0068 A' 2 , while the other values are the same as in Fig. 5c. 

4. Summary 

Computer simulations of CNTs in an EVA matrix were performed in order to 
evaluate the effect of various parameters on thermal conductivity of the composite. In 
the atomistic-based modeling, interfacial thennal resistance values, Rk, were calculated 
as a function of several functionalization parameters. The values were found to be 
strongly influenced by the grafting density, cr, and chain length, n, of short linear 
hydrocarbon chains covalently bonded to the CNT surface. Increasing chain length and 
grafting density lowers the interfacial thermal resistance. The presence of defects 
adjacent to grafting sites causes very little change in the predicted values of Rk. 

In the analytical modeling, the TC of the composite was found to be strongly 
dependent on CNT length, CNT volume fraction, CNT thennal conductivity, and 
interfacial thermal resistance. The dependence of TC on Rk was most pronounced at 
shorter CNT lengths. CNT length was a critical factor in the resulting TC of composites 
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at finite values of Rk. The TC increased with increasing cr and n, but did not exhibit a 
linear dependence. 
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Table 1 Interfacial thermal resistance values from Figs. 4a and 4b in units of m"K/W 

o 

(xlO ). The grafting density is shown across the top row. The length of the grafted linear 
hydrocarbon is shown across the left column. The values for Fig. 4b with the defects 
simulation are shown in parentheses. 
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Grafting Density of 0.01352 A " 2 , 20 chains, (6,6) SWCNT 


Fig. 1. A single walled carbon nanotube (SWCNT) with short end-grafted chains is 
depicted. 
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Fig. 2. The decay of the temperature difference between the carbon nanotube (CNT) and 
the polymer matrix. 
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Fig. 3a. The radial density profiles for the nanocomposite system with C 12 H 25 chains 
end-grafted to the nanotube in a poly(ethylene vinyl acetate) matrix for a range of 
grafting densities. The large peak at -5 A is due to the carbon nanotube wall. At high 
grafting density (cr= 0.054 A' 2 ), the profile has a qualitative change in shape. 
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Fig. 3b. The radial density profiles for the nanocomposite system with C 12 FF 5 chains 
end-grafted to the nanotube in a poly(ethylene vinyl acetate) matrix for a range of 
grafting densities. The large peak at -5 A is due to the carbon nanotube wall. The 
density profiles are due to the end-grafted chains. The density profile contribution from 
the matrix has been removed. A qualitative change in the profile occurs at high grafting 
density (cr= 0.054 A' 2 ), as the chains become extended due to crowding. 
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Fig. 4a. The interfacial thermal resistance (Rk) is plotted as a function of the grafting 
density, cr, for several end- grafted chain lengths. 
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Grafted chain 



Fig. 4b. Defects are approximated by breaking two of the bonds adjacent to a grafting 
point. 
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Grafting Density (A' 2 ) 


Fig. 4c. The interfacial thermal resistance (Rk) is plotted as a function of the grafting 
density, cr, for pristine nano tubes and nano tubes with defects. Results are shown from 
MD simulations with grafted chain length, n = 12. 
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Fig. 5a. The thermal conductivity of the nanocomposite as a function of interfacial 
thermal resistance ( R K ) for a range of volume fractions (f= 1-5%). A value of 10‘ 6 m (1 
micron) is used for the nano tube length, L. A value of 6000 W/mK is used for the 
conductivity of the nanotube, K c . The plotted points correspond to values obtained form 
the MD simulation for end-grafted chains of length n = 8. 
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Fig. 5b. The thermal conductivity of the nanocomposite as a function of interfacial 
thermal resistance ( R K ) for a range of volume fractions (/ = 1-5%). A value of 10‘ 5 m (10 
microns) is used for the nanotube length, L. A value of 6000 W/mK is used for the 
conductivity of the nanotube, K c . The plotted points correspond to values obtained form 
the MD simulation with end-grafted chains of length n = 8. 
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Fig. 5c. The thermal conductivity of the nanocomposite as a function of interfacial 
thermal resistance ( R K ) for a range of volume fractions (/ = 1-5%). A value of 10 4 m 
(100 microns) is used for the nano tube length, L. A value of 6000 W/inK is used for the 
conductivity of the nanotube, K c . The plotted points correspond to values obtained form 
the MD simulation with end-grafted chains of length n = 8. 
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Fig. 6. The thermal conductivity of the nanocomposite as a function of nanotube length 
(L) for a range of thermal interfacial resistance Rk) values. A value of 6000 W/mK is 
used for the conductivity of the nanotube, K c . The plotted points correspond to values of 
L used in Figs. 5a-c. 
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Fig. 7. The thermal conductivity of the nanocomposite as a function of interfacial 
thermal resistance ( R K ) for a range of volume fractions (/ = 1-5%). A value of 10 4 m 
(100 microns) is used for the nano tube length, L. A value of 3000 W/inK is used for the 
conductivity of the nanotube, K c . The plotted points correspond to values obtained from 
the MD simulation with end-grafted chains of length n = 8. 
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Fig. 8. The thermal conductivity of the composite to the matrix as a function of grafting 
density (cr) for a range of volume fractions (/= 1-5%). A value of 10" 4 m (100 microns) 
is used for the nanotube length, L. A value of 6000 W/mK is used for the conductivity of 
the nanotube, K c . The plotted points correspond to values obtained from the MD 
simulation with end-grafted chains of length n = 8. 
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Fig. 9. The thermal conductivity of the nanocomposite as a function of length of the 
grafted chain for a range of volume fractions (f= 1-5%). A value of 10' 4 m (100 microns) 
is used for the nanotube length, L. A value of 6000 W/mK is used for the conductivity of 
the nanotube, K c . 
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